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Measurement  of  the  Unfrozen  Water  Content  of  Soils 

Comparison  of  NMR  and  TDR  Methods 


MICHAEL  W.  SMITH  AND  ALLEN  R.  TICE 


INTRODUCTION 

The  two  methods  that  are  the  subject  of  this 
paper— pulsed  nuclear  magnetic  resonance  (NMR) 
and  time  domain  reflectometry  (TDR)— offer 
complementary  attributes  in  the  measurement  of  - 
soil  water  content.  NMR  provides  a  fast  and  accu¬ 
rate  method  for  use  with  laboratory  samples, 
while  TDR  can  be  used  for  both  in-situ  measure¬ 
ments  in  the  field  and  laboratory  experiments. 

Time  domain  reflectometry  (TDR)  measures  the 
travel  time  of  a  megahertz  pulse  through  the  soil, 
from  which  an  apparent  dielectric  constant  is  de¬ 
termined.  The  unfrozen  water  content  is  inferred 
from  a  calibration  curve  of  the  apparent  dielectric 
constant  vs  volumetric  water  content,  determined 
by  experiment.  Previously,  precise  calibration  of 
the  technique  for  determining  water  content  was 
hindered  by  the  lack  of  a  reference  comparison 
method,  which  NMR  now  offers.  The  calibration 
procedure  involves  simultaneous  measurements  of 
unfrozen  water  content  (NMR)  and  dielectric  con¬ 
stant  (TDR)  for  a  variety  of  soil  specimens.  Before 
describing  the  experiment,  we  will  briefly  review 
the  two  methods. 

Nuclear  magnetic  resonance  (NMR) 

All  atomic  nuclei  have  magnetic  moments  but 
the  quantum  energy  levels  are  characteristic  of  the 
particular  nuclear  species.  When  radio  frequencies 
are  applied,  atoms  absorb  a  certain  amount  of 
energy  to  realign  to  another  stable  position  within 
the  magnetic  field.  If  a  soil-water  mixture  is  placed 
in  a  pulsed  NMR  analyzer  and  a  single  radio  fre¬ 
quency  pulse  is  applied,  a  voltage  (which  corre¬ 
sponds  to  the  number  of  atoms  absorbing  energy) 
is  induced  in  a  receiver  coil  that  surrounds  the 
specimen.  This  voltage  is  detected  by  the  NMR 
analyzer;  its  magnitude  (minus  the  background)  is 
directly  proportional  to  the  amount  of  water  (hy¬ 
drogen)  in  the  mixture.  Therefore,  the  NMR  can 


be  used  as  a  soil  water  detector  (subject  to  certain 
considerations,  e.g.  see  Tice  and  Oliphant  1984). 

A  drop  in  signal  intensity  is  observed  as  soil 
water  freezes,  since  the  NMR  used  (Praxis  PR-103) 
was  tuned  to  the  hydrogen  proton  associated  with 
liquid  water  (10.72  MHz).  The  signal  associated 
with  the  protons  of  the  solid  ice  and  soil  constitu¬ 
ents  is  not  recorded.  Tice  et  a*.  (1982)  describe  the 
NMR  technique  for  determining  the  unfrozen 
water  content  of  soils,  and  demonstrate  the  accu¬ 
racy  of  the  technique  by  the  agreement  of  NMR 
results  with  physical  desorption  data. 

Time  domain  reflectometry  (TDR) 

TDR  is  a  common  method  for  obtaining  fre¬ 
quency-dependent  permittivity  (dielectric  con¬ 
stant)  values  in  the  VHF  through  microwave  fre¬ 
quency  range,  utilizing  Fourier  transformations  of 
picosecond-duration  pulses  incident  and  reflected 
from  a  test  material  contained  in  a  suitable  wave¬ 
guide  (e.g.  Delaney  and  Arcone  1984).  In  a  differ¬ 
ent  application,  Topp  et  al.  (1980)  proposed  ob¬ 
taining  single  (frequency-independent)  permittiv¬ 
ity  values  for  soils  using  a  non-Fourier  approach, 
in  which  the  travel  time  of  an  electromagnetic 
pulse,  launched  along  a  coaxial  line  of  known 
length  L  containing  a  soil  specimen,  is  measured. 
(Alternatively,  a  pair  of  parallel  metallic  rods  may 
be  embedded  in  the  soil  specimen.)  The  reflections 
originating  at  the  front  and  back  (beginning  and 
end)  of  a  sample  (or  in-situ  sample  volume)  are  re¬ 
corded  and  the  time  delay  t  between  these  reflec¬ 
tions  is  measured  (Fig.  1). 

Electromagnetic  theory  shows  that  the  velocity 
of  propagation  v  can  be  expressed  in  terms  of  an 
apparent  dielectric  constant  Kt  for  the  test  medi¬ 
um  and  the  velocity  of  electromagnetic  radiation 
in  free  space,  c  (3  x  10*  m  s'1): 

=  c/Jka. 


V 


(1) 


@U  =  /(*  a)- 


Figure  I.  Typical  TDR  traces  at  a  variety  of  freez¬ 
ing  temperatures. 


Since 

v  =  L/t  (2) 

then 

Aa  =  (ct/L)\  (3) 

However,  since  wet  soils  ate  generally  dispersive  at 
frequencies  in  the  low  gigahertz  range  (Hoekstra 
and  Delaney  1974,  Delaney  and  Arcone  1984),  the 
value  of  Aa  determined  by  the  travel  time  method 
will  approach  the  true  value  of  A  only  if  the  TDR 
pulse  contains  no  such  high  frequencies  by  the 
time  it  returns  from  the  end  of  the  transmission 
line  (see  also  Arcone  and  Wills  1986).  Two  fea¬ 
tures  of  the  TDR  measurement  system  used  in  our 
experiments  contribute  to  this  condition: 

1. The  TDR  model  used  (Tektronix  1502)  has  a 
relatively  slow  rise  time  (about  140  ps),  and 
thus  the  incident  pulse  has  virtually  no  fre¬ 
quency  content  above  about  800  MHz. 

2.  When  transmission  lines  of  about  20  cm  or 
larger  are  used,  most  of  the  high  frequency 
content  of  the  signal  is  dispersed  along  the 
line  and  the  reflection  at  the  end  of  the  line  is 
dominated  by  the  lower  frequencies  in  the 
launched  signal. 


EXPERIMENTAL  RATIONALE 

Topp  et  al.  (1980)  determined  an  empirical  rela¬ 
tionship  between  A,  and  the  volumetric  water  con¬ 
tent  (0u)  that  is  nearly  independent  of  soil  type, 
density,  temperature  and  salinity.  As  the  water 
content  increases,  the  time  delay  t  increases  (i.e. 
A,  increases);  thus  a  relationship  exists  as  follows: 


Patterson  and  Smith  (!<J81)  argued  that  Topp’s  re¬ 
lationship  could  be  applied  as  a  first  approxima¬ 
tion  to  the  determination  of  the  volumetric  unfro¬ 
zen  water  content  0U  of  frozen  soils,  since  the  low 
frequency  dielectric  constant  of  ice  (3.2)  is  close  to 
the  value  for  soil  minerals  (3  to  4).  This  hypothesis 
was  confirmed  by  combined  TDR  and  dilatometer 
measurements  (Patterson  and  Smith  1981,  Smith 
and  Patterson  1984),  although  results  were  re¬ 
stricted  to  temperatures  above  -3°  or  -4°C,  be¬ 
cause  dilatometer  errors  are  cumulative  with  cool¬ 
ing.  These  authors  expected,  however,  that  the  re¬ 
lationship  for  frozen  soils  would  depart  from  the 
Topp  curve  at  lower  temperatures  (lower  unfrozen 
water  contents  and  higher  ice  contents),  as  ice  is 
“substituted”  for  air  (A  =  1)  in  the  soil.  In  addi¬ 
tion,  Topp  et  al.  (1980)  passed  their  calibration 
curve  through  a  value  of  80  for  bulk  water,  al¬ 
though  their  maximum  soil  water  contents  were 
below  about  50%.  Our  measurement  system  yield¬ 
ed  a  value  for  the  apparent  dielectric  constant  of 
water  that  was  consistently  lower  than  80  (about 
72).  Thus  a  departure  from  the  Topp  curve  is  also 
to  be  expected  at  high  soil  water  contents. 

These  considerations  pointed  to  the  need  for  a 
calibration  dedicated  to  frozen  soils,  and  the 
NMR  technique  provides  a  suitable  reference  meth¬ 
od  for  this.  The  TDR/NMR  approach  offers  con¬ 
siderable  scope  for  comparison,  accommodating  a 
wide  range  of  soil  types  and  temperature  (unfro¬ 
zen  water  content). 

Further  details  concerning  the  NMR  and  TDR 
methods  can  be  found  in  Tice  et  al.  (1978,  1982, 
1984),  Topp  et  al.  (1980)  and  Patterson  and  Smith 
(1981). 


EXPERIMENTAL  METHODS 

Separate  specimens,  from  the  same  prepared 
soil  sample,  were  subjected  to  NMR  and  TDR 
measurements.  To  facilitate  filling  of  the  TDR  co¬ 
axial  specimen  container,  all  the  samples  were  sat¬ 
urated  with  distilled  water  to  the  point  of  dila- 
tancy.  The  containers  were  from  12  to  20  cm  long, 
and  1 .35  cm  in  diameter.  The  specimen  volume  was 
from  16  to  25  cm’,  and  the  specimen  mass  was  typ¬ 
ically  25  to  40  g.  The  NMR  specimens  were  com¬ 
pacted  in  test  tubes  to  the  same  density  as  in  the 
TDR  tests;  the  specimen  was  typically  12  to  15  g. 

The  two  sets  of  specimens  were  placed  in  preci¬ 
sion  temperature  baths  containing  an  ethylene 
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Table  1.  Values  of  specific  surface  area,  gravimetric 
and  volumetric  water  contents  and  dry  density  for 
soils  used  in  this  study.  Line  separates  original  17  soils 
from  those  used  to  confirm  the  Kt  vs  $a  relationship. 


Soil 

•SS/t 

w 

e<t 

West  Lebanon  gravel 

15 

0.240 

0.375 

1.563 

Castor  sandy  loam 

— 

0.261 

0.385 

1.475 

Manchester  silt 

18 

0.303 

0.432 

1.425 

Kaolinite  (KGa-1) 

23 

0.396 

0.500 

1.262 

Chena  Hot  Springs  silt 

40 

0.284 

0.414 

1.456 

Leda  clay 

58 

0.333 

0.470 

1.412 

Mayo  silty  clay 

— 

0.378 

0.495 

1.311 

Morin  clay 

60 

0.362 

0.472 

1.305 

O’Brien  clay 

61 

0.381 

0.515 

1.352 

Goodrich  clay 

68 

0.362 

0.467 

1.289 

Tuto  clay 

78 

0.733 

0.674 

0.920 

Sweden  VFB  478  clay 

113 

0.491 

0.548 

1.116 

Suffield  silty  clay 

148 

0.340 

0.455 

1  339 

Frederick  clay 

159 

0.426 

0.492 

1.209 

Ellsworth  clay 

184 

0.372 

0.450 

1.2(0 

Regina  clay 

291 

0.595 

0.572 

0.961 

Umiat  bentonite 

714 

2.25 

0.821 

0.365 

Lanzhou  silt  (saline,  China) 

34 

0.220 

0.364 

1.655 

Niagara  silt 

37 

0.222 

3.365 

1.645 

Norway  LE-1  clay 

52 

0.337 

0.477 

1.415 

Kaolinite  no.  7 

72 

0.570 

0.587 

1.029 

Athena  silt  loam 

83 

0.352 

0.456 

1.296 

Sweden  CTH  201  clay 

106 

0.508 

0.573 

1.128 

Hectorite 

419 

1.023 

0.722 

0.706 

Volcanic  ash 

474 

0.668 

0.602 

0.901 

glycol-water  mixture.  They  were  cooled  to  be¬ 
tween  -10°  and  -15°C  and  progressively  warmed 
to  0°C.  NMR  and  TDR  readings  were  taken  at 
predetermined  temperatures  on  the  warming  cy¬ 
cle.  Readings  were  also  taken  for  the  completely 
thawed  specimens.  The  TDR  traces  were  recorded 
on  a  Hewlett-Packard  7045B  X-Y  recorder;  this 
allowed  the  travel  time  to  be  determined  to  within 
±  lVo. 

Anderson  and  Tice  (1972)  demonstrated  that 
the  unfrozen  water  content  of  frozen  soils  depends 
not  only  on  the  temperature  but  also  the  specific 
surface  area  (SSA,  mJ  g  ')  of  the  soil.  Thus  for  the 
present  experiments,  a  variety  of  soils  (17  in  num¬ 
ber)  were  selected,  which  covered  a  representative 
range  in  specific  surface  area  (Table  1).  An  addi¬ 
tional  eight  soils  (Table  1),  which  were  not  part  of 
the  calibration  set,  were  also  tested  as  a  means  of 
verifying  the  relationship  determined  for  Kt  vs  0U. 


RESULTS  AND  DISCUSSION 

With  the  dielectric  constant  at  each  temperature 
determined  by  the  pulse  reflection  method  (TDR) 


Temperature  (°C) 


Figure  2.  Water  content,  NMR  and  TDR  rela¬ 
tionships  for  Leda  clay. 


and  the  unfrozen  water  content  determined  by  the 
NMR  technique,  a  typical  set  of  results  is  shown  in 
Figure  2.  Correlation  of  the  dielectric  constant 
with  the  amount  of  unfrozen  water  present  was 
achieved  by  matching  the  two  curves  on  the  basis 
of  temperature.  Since  the  NMR  values  are  on  a 
gravimetric  basis,  the  volumetric  unfrozen  water 
content  in  the  TDR  specimen  at  any  temperature 
was  obtained  by  multiplying  the  NMR  value  by 
the  dry  density  Qd  of  the  TDR  specimen: 
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Figure  3.  Water  content  vs  Ka  (data  for  all 
soils). 


e  u 


Figure  4.  Water  content  vs  K0  (data  for  three 
soils). 


6U  =  ed 

where  ew  is  the  density  of  water.  The  gravimetric 
water  content  and  soil  density  were  determined  at 
the  completion  of  each  test. 

The  comparison  procedure  resulted  in  a  plot  of 
vs  da  (Fig.  2c),  providing  a  calibration  of  the 
TDR  technique  for  unfrozen  water  content  deter¬ 
mination.  In  this  case,  a  cubic  polynomial  fits  the 
data  with  a  standard  error  (SE)  of  0.5%  in  0U, 
which  is  indicative  of  the  accuracy  possible  with 
the  TDR  technique  when  it  is  calibrated  for  a 
single  soil.  However,  the  real  purpose  of  the  ex¬ 
periment  was  to  examine  the  dependence  or  other¬ 
wise  of  the  TDR  technique  on  soil  type. 

To  this  end,  the  comparison  procedure  was  re¬ 
peated  for  all  17  test  soils,  yielding  over  220  points 
for  Kt  vs  0a  (Fig.  3).  A  third-degree  polynomial 
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Water  Content  (NMR) 

Figure  5.  Results  of  verification  tests  of  eq  4. 


fitted  to  these  data  gives  the  following  general  re¬ 
lationship: 

du  =  -  1 .458£  -  1  +  3 .868£  -  2Ka 

-8.502£-4AaJ  +  9.920£-6/f1j  (4) 

with  a  standard  error  of  estimate  of  ±1.55%  in 
0U.  The  main  reason  for  the  scatter  in  Figure  3  is  a 
subtle  but  apparent  dependency  on  the  texture 
(specific  surface  area?)  of  the  soil.  This  is  demon¬ 
strated  in  Figure  4,  where  data  for  three  soils  are 
illustrated;  it  is  evident  that,  for  a  given  water  con¬ 
tent,  the  apparent  dielectric  constant  decreases 
with  the  fineness  of  the  soil.  It  is  suggested  that, 
since  the  dielectic  constant  of  adsorbed  water  must 
be  (much)  lower  than  that  of  bulk  water,  the  large 
proportion  of  adsorbed  water  in  materials  with 
high  SSA  lowers  the  overall  dielectric  constant 
(travel  time)  of  the  soil  system  at  a  given  water 
content. 

For  most  practical  purposes,  however,  eq  4, 
which  does  not  distinguish  between  soil  types,  may 
be  used  to  determine  0U  from  the  measurement  of 
Kt  to  within  ±0.03  m’ m'*  (2  SEs),  with  a  95% 
level  of  confidence.  As  a  verification  of  this  rela¬ 
tionship,  a  further  eight  soils,  which  did  not  form 
part  of  the  calibration  set,  were  tested.  These  re¬ 
sults  form  an  independent  assessment  of  eq  4  for 
determining  the  unfrozen  water  content  of  soils. 
The  comparison,  which  is  shown  in  Figure  5,  indi¬ 
cates  that  the  general  relationship  is  satisfactory, 
and  that  the  TDR  method  can  be  used  without  re¬ 
gard  to  soil  type.  The  slight  bias  noticeable  in  Fig¬ 
ure  5  is  probably  due  to  the  phenomenon  illustrat¬ 
ed  in  Figure  4. 


Figure  6.  Comparison  of  NMR  and  TDR  data 
for  hectorite. 


The  most  notable  departure  from  the  general  re¬ 
lationship  described  in  eq  4  was  observed  in  the  re¬ 
sults  for  two  colloidal  soils — hectorite  (SSA  =  419) 
and  volcanic  ash  (SSA  =  474).  In  both  cases,  the 
unfrozen  water  contents  determined  with  eq4  were 
consistently  lower  than  the  NMR  data,  except  very 
near  0°C  (Fig.  6).  These  soils  presumably  repre¬ 
sent  an  extreme  case  of  the  phenomenon  discussed 
above.  Further  investigations  of  this  are  warrant¬ 
ed  since  it  now  appears  that  a  separate  calibration 
would  be  necessary  when  using  the  TDR  with  such 
materials. 


ERRORS 

In  light  of  the  comments  in  the  introductory  sec¬ 
tion  of  the  paper  about  KM  and  K,  it  appears  that 
use  of  a  measurement  system  distinctly  different 
from  that  used  in  the  calibration  experiments— 
such  as  TDR  equipment  with  faster  rise  times,  or 
the  use  of  very  short  transmission  lines— could 
lead  to  a  larger  error  than  suggested  by  eq  4. 

Finally,  one  other  possible  source  of  error  is  de¬ 
rived  from  observed  variations  in  the  time-base 
calibration  of  different  1502  TDR  units.  Using 
different  machines,  ve  have  measured  values 
for  distilled  water  ranging  from  69  to  77.  The  unit 
that  was  used  in  determining  eq  4  gave  a  value  of 
72  for  water  at  20°C,  in  a  20-cm  coaxial  line  (500  0 
in  air).  If  one  were  to  obtain  a  value  of  75,  say,  on 
another  unit,  then  all  measurements  using  that 
unit  should  be  reduced  by  4W  in  order  to  use  eq  4 
without  additional  error. 


CONCLUSIONS 

It  appears  that  the  TDR  technique,  which  is  ap¬ 
propriate  for  use  in  undisturbed  soil  conditions  (in 
the  field  or  laboratory),  does  not  require  special 
calibration  for  most  practical  p«n  ses.  Overall, 
the  results  obtained  from  the  cot  aed  NMR  and 
TDR  experiments,  based  on  25  soils  covering  a 
wide  range  of  SSAs  likely  to  be  normally  encoun¬ 
tered  in  the  field,  yield  a  relationship  between  ap¬ 
parent  dielectric  constant  and  unfrozen  water  con¬ 
tent  that  is  only  very  slightly  dependent  on  soil 
type.  Use  of  this  relationship  implies  a  likely  error 
of  about  ±0.03  cm5  cm'’  in  water  content  with  the 
TDR  method.  The  only  notable  exceptions  ob¬ 
served  so  far  are  the  results  for  two  soils  with  very 
high  SSAs,  for  which  the  predicted  water  contents 
were  markedly  lower  than  the  corresponding 
NMR  values. 
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